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The contribution of ground state correlations (GSC) to the non-mesonic weak decay of and 
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other medium to heavy hypernuclei is studied within a nuclear matter formalism implemented 
in a local density approximation. We adopt a weak transition potential including the exchange 
of the complete octets of pseudoscalar and vector mesons as well as a residual strong interaction 
modeled on the Bonn potential. Leading GSC contributions, at first order in the residual strong 

<N 

interaction, are introduced on the same footing for all isospin channels of one- and two-nucleon 

ON " 

induced decays. Together with fermion antisymmetrization, GSC turn out to be important for an 
accurate determination of the decay widths. Besides opening the two-nucleon stimulated decay 
channels, for j^ 2 C GSC are responsible for 14% of the rate T\ while increasing the T n /T p ratio by 
4%. Our final results for fC are: T NM = 0.98, T n /T p = 0.34 and r 2 /T NM = 0.26. The saturation 
property of Tnm with increasing hypernuclear mass number is clearly observed. The agreement 
with data of our predictions for Tnmj r n /Tp and T 2 is rather good. 

PACS numbers: 21.80.+a, 25.80.Pw. 
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I. INTRODUCTION 



The study of nuclear systems with strangeness is a relevant question in modern nuclear 
and hadronic physics [1], which also implies important links with astrophysical processes 
and observables as well as with QCD, the underlying theory of strong interactions. Vari- 
ous strange nuclear systems can be studied in the laboratory, ranging from hypernuclei and 
kaonic nuclei to exotic hadronic states such as strangelets, if-dibaryons and pentaquark 
baryons. Strangeness production can also be investigated in relativistic heavy-ion collision 
experiments, whose main aim is to establish the existence of a quark-gluon plasma. More- 
over, the cold and dense matter contained in neutron stars is expected to be composed 
by strange hadrons, in the form of hyperons and Bose-Einstein condensates of kaons, and 
eventually by strange quark matter for sufficiently dense systems. 

The existence of hypernuclei — bound systems of non-strange and strange baryons — 
opens up the possibility to study the hyperon-nucleon and hyperon-hyperon interactions 
in both the strong and weak sectors. In turn, such interactions are important inputs, for 
instance, when investigating the macroscopic properties (masses and radii) of neutron stars. 
The best studied hypernuclei contain a single A-hyperon. In a nucleus the A can decay 
by emitting a nucleon and a pion (mesonic mode) as it happens in free space, but its 
(weak) interaction with the nucleons opens new channels which are indicated as non-mesonic 
decay modes (for recent reviews see Refs. [2-6]). These are the dominant decay channels of 
medium-heavy nuclei, where, on the contrary, the mesonic decay is disfavoured by the Pauli 
blocking effect on the outgoing nucleon. In particular, one can distinguish between one- and 
two-body induced decays, AN — > nN and ANN — > nNN. The hypernuclear lifetime is given 
in terms of the mesonic (I'm = IV- + T^o) and non-mesonic decay widths (IT^m = Ti + ^2) 
by t = H/Tt = h/ \Vyi + Tnm] ■ The various isospin channels contribute to the one- and two- 
nucleon induced non-mesonic rates as follows: I\ = T n + T p = T(An — > nn) + T(Ap — > np) 
and T 2 = T nn + T np + T pp = T(Ann — > nnn) + T(Anp — > nnp) + T(App — > npp). 

One should note that, strictly speaking, the only observables in hypernuclear weak decay 
are the lifetime r, the mesonic rates r^- and T w a and the spectra of the emitted particles 
(nucleons, pions and photons). None of the above non-mesonic partial decay rates (r„, T p , 
T np , etc) is an observable from a quantum-mechanical point of view. Each one of the possible 
elementary non-mesonic decays occurs in the nuclear environment, thus subsequent final 
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state interactions (FSI) modify the quantum numbers of the weak decay nucleons and new, 
secondary nucleons are emitted as well: this prevents the measurement of any of the non- 
mesonic partial decay rates. Instead, the total width T T can be measured: being an inclusive 
quantity, for such a measurement one has to detect any of the possible products of either 
mesonic or non-mesonic decays (typically protons from non-mesonic decays). The fact that 
the detected particles undergo FSI does not appreciably alters the lifetime measurement, 
since strong interactions proceeds on a much shorter time scale than weak decays, and 

T measured = r _|_ r strong ^ r = fr/T^. 

In order to achieve a proper knowledge of the various decay mechanisms (in particular 
of the strangeness-changing baryon-baryon interactions), a meaningful comparison between 
theory and experiment must be possible. The above discussion shows that such a comparison 
requires the introduction of non-standard theoretical definitions for the non-mesonic par- 
tial decay rates (which, as mentioned, are not quantum-mechanical observables) together 
with the corresponding experimental methods for determining these rates. In our opin- 
ion, this point has not been adequately addressed in previous works and, among others, it 
has impacted on the well-known puzzle on the ratio T n /T p between the neutron- and the 
proton-induced non-mesonic rates. 

In order to explain how the total non-mesonic rate can be determined in an experi- 
ment, we have to discuss first the measurement of the mesonic rates. The pion and nucleon 
emitted in a mesonic decay both have a momentum of about lOOMeV/c. Nucleons of a 
few MeV kinetic energy cannot be observed as they are below the experimental detection 
thresholds. Mesonic decays are thus identified by measuring pions (7r~'s or it — > 77 de- 
cays). The mesonic width r„— (r^-o) is determined from the observed n~ (tt° — > 77) energy 
spectra and the total width Tt- For instance, T e * p = (A^ 7r -/A^h yp )r^ cp , N^- being the total 
number of detected 7r~'s and iVhy P the total number of produced hypernuclei. Both these 
numbers are corrected for the detection efficiencies and the detector acceptances implied 
in the measurements. The mesonic rates measured in this way thus include the effect of 
in-medium pion renormalization. Theoretical models [7, 8] also taking into account dis- 
torted pion waves obtained mesonic widths in agreement with the experimental values (in 
particular, the importance of the pion wave-function distortion was first demonstrated in 
the works of Ref. [7]). 

The experimental total non-mesonic rate is then obtained as the difference between the 
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total and the mesonic rates, = r^ p — r^ p - The experimental determination of T n /T p is 
much more involved. Indeed, this ratio must be extracted from the nucleon emission spectra, 
and this requires some theoretical input [10, 11]. FSI are very important for the non-mesonic 
processes and nucleons which have or have not suffered FSI are indistinguishable between 
each other. A theoretical simulation of nucleon FSI is thus needed and, in principle, a 
coherent sum of both kinds of nucleons must be considered when evaluating the spectra. 
Generally, FSI are accounted for by an intranuclear cascade model [9], which is a semi- 
classical scheme. 

In the present work we study the non-mesonic weak decay of hypernuclei ranging from j^B 
to ^ 08 Pb by using a nuclear matter approach implemented in a local density approximation. 
All the possible isospin channels for one- and two-body induced mechanisms are included in 
a microscopic approach based on the evaluation of Goldstone diagrams. The partial decay 
rates are derived by starting from a two-body weak transition potential. In particular, we 
investigate the effect of ground state correlations (GSC), i.e., the contribution of nucleon- 
nucleon correlations in the hypernucleus ground state. Leading order GSC contributions 
will be introduced on the same ground for one- and two-nucleon induced processes for the 
first time. The general formalism we adopt was established in Refs. [12, 13]. The weak 
transition potential for the nucleon-nucleon strong interaction contributing to the GSC we 
adopt a Bonn potential with the exchange of 7r, p, a and oj mesons. 

The paper is organized as follows. In Section II we start with general considerations 
about FSI, the definitions we employ for the weak decay rates as well as the method usually 
employed for the determination of T n /T p from data on nucleon spectra. In Section III we 
present and discuss the general framework for the evaluation of the one- and two-nucleon 
induced decay widths with the inclusion of GSC. In Section IV we make some further 
considerations about the evaluation of the widths and we discuss some former work on the 
subject. Explicit expressions for the considered GSC diagrams contributing to the one- 
nucleon induced rates are given in Section V and in Appendix A. Then, in Section VI we 
present our results and finally in Section VII some conclusions are given. 
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II. PRELIMINARY CONSIDERATIONS ON FSI EFFECTS AND ON THE DE- 
TERMINATION OF THE WEAK DECAY RATES 



The r„/T p ratio is defined as the ratio between the total number of primary (i.e., weak 
decay) neutron-neutron and neutron-proton pairs, N™^ and A^ d , emerging from the pro- 
cesses An — > nn and Ap — > np, respectively. Due to nucleon final state interactions and 
two-body induced decays, the following inequality is expected for the observables nn and 
np coincidence numbers, N nn and N np [10] [35]: 

F — Mwd I AT V ) 

P ly np JV "P 

Only N nn /N np is a quantum-mechanical observable: generally, its measurement is affected by 
thresholds on the nucleon energy and the pair opening angle [6, 14, 15]. Theoretical models 
are thus required to determine the "experimental" value of T n /T p from a measurement of 
Nnn/N n p. This unusual procedure to determine (r„/r p ) cxp makes complete sense provided 
different models are at disposal and lead to the same extracted ratio: only in such a case 
one is allowed to define this value as the experimental result for T n /T p . It is thus important 
to explore the predictions of alternative models when applied to the analysis of data. In 
the present section we go deeper into questions of this kind to show some ambiguities which 
need to be emphasized for a meaningful comparison between theory and experiment. 

Let us first illustrate in some detail the procedure normally adopted to extract (T n /T p ) cxp 
from measurements of N nn /N np [10, 11]. Each one of the non-mesonic weak decay channel 
takes place by the emission of two or three primary nucleons. These nucleons propagate 
within the nuclear environment and cannot be measured. The strong interactions with 
the surrounding nucleons can change the charge and the energy-momentum of the primary 
nucleons; some of them can be absorbed by the medium and the emission of additional 
(secondary) nucleons can occur as well. All these processes are generically designated as 
final state interactions (FSI): they do not have to be included when calculating the decay 
rates, but the observable nucleon spectra, i.e., N nn and N np , are crucially affected by them. 
One has to emphasize that, on the contrary, baryon-baryon short range correlations in both 
the initial and the final states as well as mean field effects on the single particle wave- 
functions are genuine contributions to the decay rates. 

FSI pertain to the same quantum-mechanical problem which starts with the A decay and 
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ends with the detection of the particles emitted by the hypernucleus. In a strict quantum- 
mechanical scheme, FSI cannot thus be disentangled from the weak interaction part of the 
problem: this is an analogous way of expressing the fact that the weak decay rates are 
not measurable. However, up to now FSI have been simulated by means of semi-classical 
models, i.e., by intranuclear cascade codes (INC) [9] acting after the weak decay, thus losing 
quantum-mechanical coherence. In such INC analyses, both one- and two-nucleon induced 
decays are included as inputs and one proceeds to fit N nn /N np data in order to determine 
the value of (r„/r p ) cxp . Technically, this is achieved by applying Eq. (16) of Ref. [11] (see 
also Eqs. (1) and (2) of Ref. [10] and Eq. (3.14) of Ref. [16]), which is an exact relation 
only neglecting quantum coherence among the final, observable nucleons. Note that such a 
procedure also requires a theoretical estimate for the ratio Y 2 /Ti. In other words, present 
nucleon-nucleon coincidence data only allows us to determine a correlation property between 

r n /r p and ryrv 

In general terms, one could wonder if it is possible to identify those quantum-mechanical 
contributions whose classical limit leads to a factorization between the weak decay process 
and the INC rescattering. This is a relevant question since in the theoretical evaluation 
of the non-mesonic decay rates FSI contributions must not be included; one indeed aims 
to extract the contribution of the elementary AiV — > nN and Ann — > nNN processes by 
studying hypernuclear decay. Unfortunately, the above question does not seem to have a 
simple solution. Although we make here some considerations about this point, we believe 
that a complete answer to it goes beyond the present contribution. 

Let us illustrate, by using an example, the nature of the problem. Consider the A self- 
energy diagram (a) of Fig. 1. This is a (time-ordered) Goldstone diagram where the weak 
transition potential y AN ^ NN , which is a two-body operator, produces an intermediate 2plh 
configuration; afterwards, the action of the nucleon-nucleon strong interaction V NN creates 
a further Iplh pair and leads to a 3p2h final state. In terms of amplitudes, y AN -+ NN 
produces two nucleons, one of which then strongly interacts with another nucleon, ending 
in the emission of three nucleons. Since the potential V NN acts after v AN ^ NN , diagram (a) 
contains a FSI effect and we argue that it must not be included when evaluating the non- 
mesonic decay rate. Note that the idea of an interaction taking place after or before another 
one is a valid statement here as we are working with Goldstone diagrams. On the contrary, 
diagram (b) of Fig. 1 represents a ground state correlation (GSC) effect. It corresponds to 



6 




FIG. 1: Goldstone diagrams for FSI (a) and 2p2h GSC contributions (b) for three nucleon emis- 
sion. The dashed and wavy lines stand for the potentials y AN ^ NN anc [ yNN^ reS p ec ti V ely. The 
diagram (a) has poles on the 2p\h and 3p2h configurations, while (b) has a single pole on the 3p2h 
configuration. For the present discussion we only consider the 3p2h poles indicated by the dotted 
lines. 

an amplitude in which the A decays by interacting with a correlated nucleon pair. Since 
the nucleon-nucleon interaction takes place before the action of the weak transition, this 
diagram must be considered when evaluating the decay rate r 2 . 

Note also that the Goldstone diagrams (a) and (b) are two different time orderings of the 
same Feynman diagram. If T 2 were an observable, it would have to be evaluated by means 
of Feynman rather than Goldstone diagrams; both diagrams (a) and (6) would contribute 
to T 2 . These diagrams must actually be taken into account when evaluating the observable 
nucleon spectra. However, here we argue that, since T 2 is not an observable, some of the 
Goldstone diagram should not be included in the theoretical definition of this rate. The class 
of diagrams that does not contribute to T 2 depends on the definition one adopts for FSI. 
Our definition leaves aside those Goldstone diagrams, like diagram (a) in Fig. 1, in which 
at least one nucleon-nucleon interaction takes place after the weak transition potential. If 
on the other hand one were to include diagram (6) in the calculation of the widths, then it 
would not be clear how to identify the diagrams incorporating FSI effects. 

A similar analysis to the previous one holds for the one-nucleon induced rates. Summariz- 
ing, we assume that one- and two-nucleon induced decay widths, which are not observables, 
are interpreted in terms of Goldstone diagrams in which no FSI effect is present. All the 
Goldstone diagrams in which at least one nucleon-nucleon interaction takes place after the 
weak transition potential must not be included when evaluating the decay rates. Any Gold- 
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stone diagram representing a GSC is instead a genuine contribution to the rates. In the 
calculation of the observable nucleon spectra, a description in terms of Feynman diagrams 
must instead be employed. 

III. MANY BODY TERMS IN THE NON MESONIC DECAY RATES 

Let us consider the one and two-body induced non-mesonic weak decay width for a 
A-hyperon with four-momentum k = (ho, k) inside infinite nuclear matter with Fermi mo- 
mentum /cp. In a schematic way, one can write: 

r 1(2) (M F ) = £ \(f\V AN ^ NN \0) kF \*8(E f -E ) , (2) 

/ 

where \0)k F and |/) are the initial hypernuclear ground state (whose energy is E ) and the 
possible 2plh or 3p2h final states, respectively. The 2p\h (3p2h) final states define I\ (r 2 ). 
The final state energy is Ef and y AN -* NN j s the two-body weak transition potential. 

The decay rates for a finite hypernucleus are obtained by the local density approxi- 
mation [17], i.e., after averaging the above partial width over the A momentum distri- 
bution in the considered hypernucleus, [^(fe)! 2 , and over the local Fermi momentum, 
&f(?") = {37r 2 p(r)/2} 1 / 3 , p{r) being the density profile of the hypernuclear core. One thus 
has: 

r 1(2) = J dk\Mk)\ 2 j dr\Mr)\ 2 r im (k,k F (r)) , (3) 

where for ip A (r), the Fourier transform of ip A (k), we adopt the lsi/ 2 harmonic oscillator 
wave-function with frequency hw (= 10.8 MeV for A 2 C) adjusted to the experimental energy 
separation between the s and p A-levels in the considered hypernucleus. The A total energy 
in Eqs. (2) and (3) is given by k = m A + k 2 /{2m A ) + V A , V A (= -10.8 MeV for \ 2 C) being 
a binding energy term. 

Since y AN ~>- NN [ s a two-body operator, the emission of two nucleons is originated either 
from the Hartree-Fock vacuum or from GSC induced by the nucleon-nucleon interaction. 
At variance, the emission of three nucleons can be only achieved when y AN ~> NN ac ts over a 
GSC. It is therefore convenient to introduce the following hypernuclear ground state wave- 
function [18]: 

low = m?) ( i ) - E w ^7' )d ; e \php' h ')\ ^ | Pa) , ( 4) 
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FIG. 2: Direct (D) and exchange (E) Goldstone diagrams for the 2p2h GSC induced by the nuclear 
residual interaction V NN . 

where | ) is the uncorrelated core ground state wave-function, i.e., the Hartree-Fock vacuum, 
while the second term in the rhs represents 2p2h correlations and contains both direct (D) 
and exchange (E) matrix elements of the nuclear residual interaction V NN . Besides, \p\} is 
the normalized state of the A, the particle and hole energies are denoted by e,i and: 

-1/2 



p,h,p',h' 



(php'h'\V NN \ ) D+E 



£p — £/i + £p> — Eh' 



(5) 



is the ground state normalization function. The particular labeling of Eqs. (4) and (5) is 
explained in Fig. 2. The explicit expression for M{k F ) is given in Ref. [19]. 
By inserting Eq. (4) into Eq. (2), for Y 1 one obtains: 
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the final states |/) being restricted to 2plh states. For T 2 one has: 

r 2 (k,k F ) = M 2 {k F )Y,t{E f -E Q ) 

f 
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£p — £h + £p> — £h> 



where the final states are given by 3p2h states. Note that all the matrix elements of V NN 
and y^-^NN appear in the antisymmetrized form. 

Let us focus now on the kind of diagrams contributing to F 1 and T 2 . This discussion 
is done in terms of transition amplitudes rather than self-energies. In Fig. 3 we report 
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(a) 




(61) (62) (63) 




FIG. 3: Transition amplitudes contributing to T±. A double-line (without arrow) represents the 
A(1232) resonance. 

some of the most representative transition amplitudes which contribute to lY All diagrams 
but (a) are originated by a GSC. Only the contribution of diagram (a) to r\ has been 
calculated microscopically up to now. The line (b) represents typical 2p2h correlations. The 
contribution (c) is a contact term involving a imNN strong vertex, while line (d) represents 
the contribution of the A (1232) resonance. It should be mentioned that there has been a 
great deal of controversy around the theoretical determination of the T n /T p ratio and the 
challenging comparison with data. In these discussions, all theoretical efforts have been 
devoted to the (a) term only; the remaining ones have simply been ignored. 

A similar analysis can be done for T 2 starting from the amplitudes of Fig. 4. Again, only 
the (a) term has been evaluated up to now in microscopic calculations [13]. The graphs 
in Figs. 3 and 4 are only representative cases. For instance, also the amplitude of Fig. 5 
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(a) (b) (c) 




FIG. 4: Transition amplitudes contributing to IV 




FIG. 5: Transition amplitude contributing to Ti and involving a strong interaction V between 
the hyperon and a lplh pair. 

should be included when calculating IY Unlike the other amplitudes of Figures. 3 and 4, 
the one in Fig. 5 involves a strong interaction V AN between the A and a lplh pair (i.e., 
a lplh GSC) and then the usual action of the weak transition potential. Apart from the 
explicit calculation, such a contribution could in principle be included in an effective way 
through the calculation of diagram (a) of Fig. 3 with a suitably chosen weak transition 
potential y AN ^ NN . However, based on the absence of isovector-meson exchange in the 
strong potential V AN , one may anticipate a small effect of this amplitude. Other amplitudes 
will provide important contributions. In the graph (a) of Fig. 4 the weak transition potential 
can also be connected to a hole line [13]. In addition, since V NN and y AN ^ NN are two- 
body operators whose matrix elements are antisymmetrized, Pauli exchange terms must be 
considered as well [19]. 

All the graphs in Figs. 3 and 4 have the same initial state, which is the hypernuclear 
ground state. The final state of the graphs in Fig. 3 (Fig. 4) is a 2plh (3p2h) state. To obtain 
the various decay width, all graphs representing transitions amplitudes with the same initial 
and final states are added and then squared. For instance, from Fig. 4 one obtains a total 
of six direct diagrams: the square of each individual amplitude plus the three interference 
terms. For the amplitudes in Fig. 3 there is a total of twenty-eight different direct terms. 
In addition, antisymmetrization considerably increases the amount of diagrams. From our 



11 



previous works it is clear to us that a full microscopic evaluation of each term is mandatory 
for several reasons. First, a raw estimation of a remarkable amount of different diagrams 
makes the final result quite uncertain. Secondly, there is no ground to evaluate differently 
the diagrams originated from Fig. 3 and those from Fig. 4: once a microscopic calculation 
is performed for the square of diagrams (a) of Figs. 3 and 4, the same should be done for 
the remaining contributions, which are all leading order GSC contributions. 

In the present work, as a further step towards the calculation of the whole set of diagrams 
relevant for the non-mesonic decay, the one-nucleon induced widths originated from the 
sum of the transition amplitudes (a) plus (61) of Fig. 3 are evaluated for the first time. 
Accordingly, the two-nucleon induced rates are instead obtained from the amplitude (a) of 
Fig. 4 by following Ref. [19]. Antisymmetrization is coherently applied to all contributions. 
Before proceeding with the formal derivation of the decay widths, in the next Section we 
first point out additional observations on the evaluation of the decay rates and on previous, 
related work. 



IV. FURTHER CONSIDERATIONS ON THE EVALUATION OF THE WEAK 
DECAY RATES 

Let us start this discussion by paying attention to the twofold effect of the nuclear residual 
interaction y NN within the matrix elements of Eq. (2). When y NN acts on the uncorre- 
cted hypernuclear ground state | ), as in Eq. (4), one has a GSC. Alternatively, V NN may 
introduce medium effects on the weak transition potential v AN ~* NN . Both effects must be 
taken into account when calculating the decay rates. In addition, V NN may modify the final 
states \ f) = \2plh) or \3p2h): for instance, acting on a \2plh) final state, it can produce a 
\3p2h) state, as in Fig. 1(a); this results in a FSI which does not contribute to Eq. (2). 

Concerning the medium effects previously mentioned, let us discuss some aspects of 
the work of Ref. [20]. Here, V NN introduces medium modifications on the mesons prop- 
agators appearing in y AN ~> NN through the direct part of the RPA (ring approxima- 
tion): schematically, in our scheme one simply has to replace y AN ~* NN with y AN ^ NN — 
y AN ^ NN 1 1 1 _ nV'^l, where the polarization propagator II contains Iplh and lAl/i con- 
tributions in Ref. [20]. Note that, since only the absolute value of the ring propagator is 
kept, the modified weak transition potential remains a real function. This approach thus 
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FIG. 6: Goldstone diagrams for FSI (a) and 2p2h GSC contributions (b) and (c) leading to two 
nucleon emission. 

represents a refinement of the weak transition potential and is consistent with Eq. (2). 

We emphasize that the mere use of diagrams when discussing the formalism developed 
in Ref. [20] or the present one could be misleading. For the approximation considered in 
Ref. [20], in Eq. (2) one has to employ the matrix element (f\V AN ~* NN \0)D+E of the modi- 
fied weak transition potential (the corresponding direct and exchange self-energy diagrams 
are shown in Fig. 2 in Ref. [20]). By making an expansion of the square of this matrix 
element in the ring series, the two terms at first order in V NN correspond to a self-energy 
contribution which matches exactly with the (a) diagram in Fig. 6, where the final state 
\f) = \2plh) corresponds either to the upper or the lower bubble. Nevertheless, the same 
diagram could also be associated to the direct part of the following product of matrix el- 
ements: (f\V NN \i) d+e (i\V AN ^ NN \0)D+E, where \i) is a 2plh intermediate configuration. 
But, since this product contains a FSI, it is not a correct contribution to the decay rates 
of Eq. (2). Antisymmetry of this product of matrix elements gives rise to a total of eight 
self-energy diagrams, which are shown in Fig. 4 in Ref. [21] and used there to calculate the 
(observable) spectra of the non-mesonic weak decay nucleons. From the analytical point of 
view, the product (f\V NN \i) d+e (i\V AN ^ NN \0) d+e is clearly different from the term at first 
order in V NN entering (/|^ A ^^ JVAr |0),D+.E. When the comparison is done using the full set 
of direct plus exchange diagrams, FSI and the medium modifications on the weak transition 
potential are manifestly different effects. Only the latter can be included in the calculation 
of the decay rates. 

As a final remark for this section, we observe that the amplitudes (a) and (61) of Fig. 3 
produce the self-energy diagrams (b) and (c) of Fig. 6. They are GSC terms and thus 
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contribute to the decay rates. Conversely, the (a) diagram of Fig. 6 must be left aside in the 
calculation, unless one considers it as a medium modification on the weak transition potential 
(but then, other medium modification contributions should be considered simultaneously), 
as done in Ref . [20] . The Goldstone diagrams of Fig. 6 are the three possible time orderings 
of the same Feynman diagram. Again, we stress that the fact that one out of three diagrams 
in Fig. 6 will not be included in our calculation of the decay rates makes sense since these 
rates are not observables and thus do not have to be described by Feynman diagrams. 



V. FORMAL DERIVATION OF THE DECAY RATES T n AND T p INCLUDING 
GSC 

In Fig. 3 we have shown a set of amplitudes which contribute to the decay rate Ti of 
Eq. (2). Only the amplitude (a) has been evaluated explicitly up to now. In the present 
work we extend the microscopic approach to include the amplitude (61), which originates 
from GSC contributions that we expect to be important. 

Before proceeding with the derivation of decay widths, it is convenient to give the expres- 
sions for the potentials. The weak transition potential y AN ^ NN and the nuclear residual 
interaction V NN read: 

v ^ NN{N N Kq)= J- TMin V%J™W(q) , (8) 

T A(JV)=°: 1 

where the isospin dependence is given by 

{ 1 f ; rTAw= ° (9) 

^ n • r 2 for r A(iV) = 1 . 

The values and 1 for ta(at) refer to the isoscalar and isovector parts of the interactions, 
respectively. The spin and momentum dependence of the weak transition potential is given 
by: 

V AN^NN {q) = {GFml) {StM (Ti ^ + sM a2 .~ + PctM (1Q) 

+Pl,t A ( ( 1) (T i qo- 2 q + P T , TA (q)(cri x O) ' (°2 x q) 
+iSv,r A (<l)( cr i x °"2) • q} , 

where the functions S TA (q), S' TA (q), P c , TA {q), Pl,t k {q), Pt,t a (q) and S v , Ta (q), which include 
short range correlations, are adjusted to reproduce any weak transition potential. 
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The corresponding expression for the nuclear residual interaction is given by: 

K N (v) = £{Vc,rM + v L ,r N (<i>i-q<T2-Q (n) 

+VT,T N {q)((Ti x q) ■ (02 x q)} , 

where the functions Vc,t n (q), VL jTJV (g) and Vt, tn (q) are also adjusted to reproduce any 
nuclear residual interaction. 

In particular, \/ A7V_>iVAr is represented by the exchange of the ir, 77, K, p, 00 and K* mesons, 
within the formulation of Ref. [22], with strong coupling constants and cut-off parameters 
deduced from the Nijmegen soft-core interaction NSC97f of Ref. [23]. For V NN we have 
used a Bonn potential [24] in the framework of the parametrization presented in Ref. [25], 
which contains the exchange of ir, p, a and u mesons. 

We give now explicit expressions for the partial decay width Ti(k, k F ) of Eq. (6), which 
for convenience is expressed in terms of its isospin components T n (k, kp) and T p (k, kp). Let 
us first rewrite Eq. (6) as follows: 

r nip) (k, kp) = r° n(p) (k, kp) + r^ sc (fc, k F ) + r«g(fc, k F ) , (12) 

where: 

r° n(p) (k,k F ) = ^(k F )J2m~E ) \(f\v^ NN \ PA ) D+F \ 2 , (is) 

r°-™ c (k,k F ) = -2AT 2 (ME E S(E f -E ){p A \(V^ NN y\f)o + E 

f p,h,p',h' 

x(flV ^N\ php > h >. PA ) D JP h P' h '> P^ VNN \Pa)d + e 



/ p,h,p',h' 



£ P — £h + e p' — e h' 

AN-^NN\ 



(f\V™^»\php'h'; p A ) D+E 



X 



(php'h'; PA \V NN \p A ) 



D+E 



£p — £h + £ p' — e h> 

The first component, r°^, is the contribution from the uncorrelated hypernuclear ground 
state, the third one, r^j, result from ground state correlations, while r°^ sc is the inter- 
ference term between correlated and uncorrelated ground states. 

It is now convenient to consider the following decomposition, dictated by the isospin 
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quantum number: 

r°n (p) (k,k F ) = r nl)(k,k F ) (14) 

P,Q=D,E 

= Y Y r T P A Q T A ,n(p) T T®T A ( k ' k F) ' 

P,Q=D,E T A /,T A =0,1 

r°<S sc (fc,M = E r^(*,*,) 

P,Q,Q'=D,E 

= Y Y 7 'r A Q TAT N ,n(p) T r^? A T N ( k , k F) , 

P,Q,Q'=D,E t a ,,t a ,t n =0,1 

P',P,Q,Q'=D,E 

^ 'T N ,T Al T A T N , Tl(p) 1 T N ,T A ,T A T N \. n 'l ) 1 

P',P,Q,Q'=D,E T N ,,T A ,,T A ,T N =0,1 

where P', P, Q, Q' = D or E refer to the direct or exchange character of the matrix elements 
of Eq. (13). The isospin factors are given by: 

r r A Q r A ,n(p) = Y ( t ^\ r A ,\f)p(f\0 TA \t A ) Q , 
f, isospin 

r r P A Q r A r N ,n(p) = Y ( t ^\ T A ,\f}p(f\ T A \tpt h t p d h ,,t A ) Q 
f, isospin 

x(t p t h t p >t h >,t A \0 TN \tA)Q> , 

^T P N P A ,? A T N ,n(p) = Y ( t ^\°r N Mpt] l tpdp l: tA)p l (tptf l tpdp l: tA\0 TAl \f)p 
/, isospin 

X (f\^T A \tpthtp'th',tA)Q{tpthtp't h i,tA\0 TN \tA)Q' , 

where the summations run over all the isospin projections t's, with the constrain that the 
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FIG. 7: Direct Goldstone diagrams corresponding to the square of the amplitude sum (a) + (61) 
of Fig. 3. See the decomposition of Eq. (14). 



emitted particles are nn for T n and np for r„. For the partial decay widths we instead find: 



M\k F ) (-1)" £ 5(£/ - £o) 
/ 

x(PA|(V^ JVAf (g')) t l/)p(/|v4™(g)bA) Q , 
-2^ 2 (M(-i) n E E 

/ p,h,p',h' 

x( P A|(v^ JV7V (g , )) t l/)p(/|v4 Ar ^ Ar (g)bV^; pa>< 



(15) 



(16) 



X 



(pV^5 PA|V^(t)|p A >g' 



V P'PQQ' (u u\ 



?p — £h + £p' — £ h' 

M\k F )(-irj2 E E s ( E f- E °) 

f p,h,p>,h> P,h,P',h' 

( PA \(V^(t')y\pXp',h'; pa)p> 



(17) 



X ■ 



x{p,h,p',h'; PA\(V^ NN (q')y\f} P 



x(/|V. 



AN—^NN 



(q)\php'h'\ p A ) ( 



x ■ 



Note that the values of the energy-momentum carried by the particles and holes lines de- 
pends on the topology of the corresponding diagram, while n is the number of crossing 
between fermionic lines. 

Let us now apply the above formalism to a model including the amplitudes (a) and (61) 
of Fig. 3. Four direct self-energy diagrams correspond to the square of the amplitude sum 
(a) + (61); they are given in Fig. 7. Note that these diagrams admits a single cut, giving rise 
to a 2plh final state. The DD diagram contributes to the partial widths Y° n ^ of Eq. (14). 
The two DDD diagrams, which have the same numerical value and are interferences between 
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FIG. 8: Goldstone diagrams for the partial rates contributing to Eq. (14). 

the amplitudes (a) and (61) of Fig. 3, are included in the partial widths r°~^j sc . Finally, 
the diagram DDDD contributes to r^^j. Many exchange diagrams are obtained from the 
antisymmetrized amplitude sum (a) + (61): one PQ exchange diagram is the partner of the 
DD one of Fig. 7; seven PQQ' exchange diagrams are companions of each one of the DDD 
ones; fifteen P'PQQ' exchange diagrams add to the DDDD one. 

Formal expressions for T° n{p) can be found in Ref. [12]. The T^^'s contributing to T°^ sc 
(see Eq. (14)) correspond to the diagrams of Fig. 8. By replacing, in Eq. (16), the sum over 
momenta by integrals and by performing the energy integrations and the spin summation, 
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the following expression for T^^ TN can be obtained: 

'4 (2tt) 8 ^ r ""*' m 



= X\k F )\\^{G F rniylt (18) 



x 



J J J dqdhdh'W^ TN (q, q ',t) 



xQ(k,q,q',t,h,ti,k F ) ]^ S(q - ( £ h'+q~ £ h'^ ' 

~ £ 2p2h 

where q = k — £fc_q — Vjv, Vn being the nucleon binding energy, while the functions 
W^®^ TN (q, q', t) and 6(/c, q, q', t, h, h', k F ) and the energy denominator are specific of 

each PQQ' contribution. The function W^^ TN (q, q', t) contains the momentum dependence 
of the nuclear residual interaction and the weak transition potentials and the spin summa- 
tion, while Q(k, q, q', t, h, h', k F ) is a product of step functions which defines the phase space 
of particles and holes. 

In the present section we present the explicit expression for the direct term T®®^ Tn ; the 
other seven ones are displayed in Appendix A. We obtain: 

?rZr N ^k F ) = M\k F ))^(G F mlY§- (19) 



x y y y ^^'< d ta u?) 

X0(q o )6(\k -q\- k F )9{\q - h\ - k F \)6(k F - \h\) 

x6(\q + ti\-k F \)6(k F - \h'\) 
1 

x- 



-*-(^-^) %0 " (£ ^" Eft ' )) - 

The expressions for Q(k, q, q', t, h, h', k F ) and are self-evident. Moreover: 



W?XtM = 8 iKAlKM + PCrJq)Pc,rM]Vc,r N (q) (20) 

+ [S T Jq)S TA (q) + P L , TA ,(q)P L , TA (q)]V L , TN (q) 

+2 [5v >TA ,(9)5v; TA (9) + PT,r Al (q)PT,rM]VT,rM} ■ (21) 

Eq. (19) can be simplified by introducing the functions: 

I(t "" q) = J^f j dh ' ^ q + h '\- kF \'> 0i - kF -\ h '\'> S{q o- E h' + q + E K) < 

KM = 1 v f ih m^hmt^m , (22) 



( 2 ^) 3 J <lQ-( £ h-q- £ h) 
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where X(g , q) is the imaginary part of the Lindhard function and the explicit expression for 
TZ(q , q) is given in Appendix B. Therefore: 



yDDD 

T A IT A T N 



x{[S' T Jq)S' TA ( q ) + Pc,rAq)Pc,rM^c,rAQ) 
+[S TAl (q)S TA (q) + P L , TA ,(g)P L , rA (g)]V L , TJV (g) 

+2 [S v , T Jq)S ViTA (q) + P T , T Jq)PT,rM}VT,rM} 

x1Z(-q ,q)l(q ,q) ■ 



(23) 



Then one has to perform the isospin summation to obtain 

Tfj ,rA,rjv=0,l 

The final result obtained after the local density approximation is therefore: 
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^DDD 



o f-pDDD , -pDDD , -pDDD , pDDDi 
z I 1 111 "r 1 000 010 101 J 



2{5r 
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i -pDDD _ -pDDD 
"T 1 000 1 010 



-r 



DDD 
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Finally, we present the partial rates corresponding to the diagram DDDD of Fig. 7. By 



applying the same procedure used for 177^ to Eq. (17) we obtain: 



^DDDD 

T M' T A' T ^ T N 



(k,k F ) = 



(2vr) 2 



■(G F ml) 2 (^) 2 / dq9(q 



) 



x6(\k -q\- k F ){(S' TA ,S' TA + Pc,t a ,Pc,M, 

+ (S TA ,S TA + Pl,t a ,Pl,t a )Vl,t n 

+2 (SV,t a ,£V,t a + Pt,t a ,Pt,t a )Vt,t n } 
x-R. 2 (-q ,q)l(q ,q) , 



TN 



(26) 



and 



DDDD 

n 

DDDD 



Af-pDDD | -pDDD | -pDDD , -pDDD\ 
'H 1 1111 "r 1 0000 ~r 1 0101 1 1010 J ) 



DDD i t^DDD t^DDD t^DDD 



= 4{5r^i +r 



0000 



0101 



1010 



(27) 



after performing the local density approximation. 

In this paper the T n ^ w exchange terms will be neglected. Indeed, from our numerical 
results discussed in the next Section it turns out that already the direct contribution r^^ DD 
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is small and approximately one order of magnitude smaller than • Moreover, according 
to the results obtained for the 's, P'PQQ' exchange contributions are expected to be 

even smaller than the direct term DDDD. 



VI. RESULTS 



In the previous Section we have seen how the neutron- and proton-induced decay widths 
can be written in the form: 

r _ pO i rO-GSC , r GSC / q\ 
l n{p) - 1 n (p) + 1 „ (p) + i n(p) 

\^ V PQ , \^ r PQQ' , V" V P'PQQ' 

P,Q=D,E P,Q,Q'=D,E P',P,Q,Q'=D,E 

-^n(p) being the rates obtained for an uncorrelated hypernuclear ground state, the rates 
originated by ground state correlations and r°~^j sc the rates resulting from the interference 
between uncorrelated and correlated ground states. 

For the present scheme containing the transition amplitudes (a) and (61) of Fig. 3, where 
antisymmetrization is considered for the weak transition potential y AN ^> NN an d the nuclear 
residual interaction V NN , we obtained: two contributions to r°^, which are = T^fp) 
and = and are generated by the square of amplitude (a); eight different 

contributions to r°~^j sc , which are interferences between the (a) and (61) amplitudes; sixteen 
different r^ > / > ^ contributions to T^^, which originate from the square of amplitude (61). 
An early evaluation of has been performed in Ref. [12], while r°~^j sc and Tnfp) are 

P 1 POO' 

discussed here for the first time. Among the r , y w 's, here we only calculate the direct 



terms T™? D . 



We start by discussing the relevance of the Pauli exchange terms in T Q n ^ and r T , . Our 
results for and are given in Table I for the decay of ]^C. Note that, for symmetry, 
■^n(p) are twice the sum of and ^nfp)- Exchange terms contribute to the uncorrelated 
rates for neutron-induced (proton-induced) decays by 5.1% (0.3%). Thus, they tend to 
increase T n /T p while having a very small effect on IY 
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TABLE I: Direct and exchange and Tp terms for j^ 2 C in units of the free A decay rate, 
r° = 2.52 • 10~ 6 eV. The first column indicates the two different isospin channels and their sum. 
Note that r D , D , = T E f , and T D f, = T EE . 

n(p) n(p) n(p) n(p) 



Channel 


2T dd 


2T de 


r° 


An — > nn 


0.146 


0.008 


0.154 


Ap — >■ np 


0.469 


0.002 


0.470 


sum 


0.615 


0.009 


0.624 



In Table II we present predictions for the T^ QQ ' and T E ®®' contributions derived from the 
Goldstone diagrams of Fig. 8, again for \ 2 C. As expected, the direct terms T EDD and Tp DD 
are the main contributions. Nevertheless, the effect of antisymmetry on the two isospin 
channels is significant: it increases r° _GSG by 34% while decreasing r° _GSC by 8%. The 
overall effect on r5~ GSC = r°~ GSC +rp" GSC is a very small increase, of 2%. We note that, with 
topologically equivalent diagrams, in Ref. [21] a similar quasi-cancellation between neutron- 
and proton-induced decays has been found in nucleon spectra calculations. Moreover, in 
Ref. [19] it has been shown that the evaluation of the GSC exchange terms is important for 
the rate T 2 as well. We emphasize that the exact evaluation of exchange diagrams has been 
mostly ignored in the literature. It is usually a quite involved (but necessary) task, given 
the rapidly increasing number of terms one has to consider when going to higher orders in 
the nuclear residual interaction. Unfortunately, there is no general rule to anticipate the 
need for the evaluation of exchange terms when the corresponding direct contribution is 
important. 

In Table III we present the different contributions to the rates T n and T p of Eq. (28). The 
uncorrelated parts T° and T° dominate over the remaining ones: Tj = T° + T° constitutes 
the 86% of the total r\. Then, r?~ GSC = ro- GSC + r°- GSC and r GSG = r GSC + r GSC represent 
13% and 1% of I\, respectively. We remind the reader that r G ^j are calculated from the 
direct diagram DDDD in Fig. 7, while P'PQQ' exchange terms are neglected. This omission 
is justified by the smallness of the direct contributions r^^ DD : the neglected exchange part 
of r GSC should contribute to Ti by less than 1%. Thus, a challenging calculation of the 
fifteen P'PQQ' exchange diagrams can be reasonably avoided. 
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TABLE II: Direct and exchange Tn®® and Tp®® terms for j^C obtained from the diagrams of 
Fig. 8. The first column indicates the two different isospin channels and their sum. 



Channel 


pDDD 


-pDDE 


YDED 


YEDD 




At?, — V Tin, 
Ap — > np 


0.022 
0.071 


-0.002 
0.005 


-0.009 
-0.027 


-0.004 
-0.011 




sum 


0.093 


0.003 


-0.036 


-0.015 




Channel 


YDEE 


YEDE 


YEED 




r 0-GSC 


An — > nn 
Ap — > np 


0.006 
-0.008 


0.008 
0.009 


0.006 
0.025 


0.002 
0.002 


0.029 
0.066 


sum 


-0.003 


0.017 


0.031 


0.004 


0.095 



TABLE III: Predictions for the one-nucleon induced decay rates of Eq. (28) for ^ 2 C. The first 
column indicates the two different isospin channels and their sum. 



Channel 


r° 


p0-GSC 


r GSC 


r 


An — > nn 
Ap — > np 


0.154 
0.470 


0.029 
0.066 


0.002 
0.008 


0.185 
0.544 


sum 


0.624 


0.095 


0.010 


0.729 



Our predictions for the one- and two-nucleon induced decay rates for ^ 2 C are given in 
Table IV and compared with the most recent data by KEK [26] and FINUDA [27]. For 
completeness, we report results without and with the inclusion of antisymmetrization and 
GSC. It should be noted that the hypernuclear ground state normalization function N(k F ) 
of Eq. (5) equally affects Ti and T 2 . This function is not identically equal to one only 
when GSC are present. Therefore, the Ti result without GSC and with exchange terms of 
Table IV, 0.74, is bigger than the prediction for of Table III, 0.62, which has been obtained 
instead by including both GSC and antisymmetrization in the normalization function. This 
comparison gives an idea of the importance of a proper normalization of the hypernuclear 
ground state. GSC produces a sizable increase in the value of T NM , thanks to the opening of 
the two-nucleon induced channel, while Ti remains practically unaffected. The effect of GSC 
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on the r„/T p ratio is a small increase of 4%, which is due entirely to the exchange terms in 
r°~ GSC and r°~ GSC (see Table II). Antisymmetrization on the contrary introduce an increase 
of r x and a reduction of T 2 , and as a result a sizable reduction oiT 2 /T 1 . We conclude that 
GSC are important to get agreement with data on Tnm, while antisymmetrization is crucial 
to reproduce the data for IT^/IY Note indeed that only with the set of results including 
both exchange terms and GSC we can achieve an overall agreement with all data. 

TABLE IV: The non-mesonic weak decay widths of ^ 2 C. Results are given without and with the 
contributions of antisymmetrization and ground state correlations. The most recent data, from 
KEK [26] and FINUDA [27], are given for comparison. 



Ant./GSC 


r n 


r P 


Ti 


r 2 






i^/Tnm 


no/no 


0.15 


0.47 


0.62 





0.62 


0.31 





yes/no 


0.18 


0.56 


0.74 





0.74 


0.33 





no/yes 


0.15 


0.47 


0.61 


0.31 


0.91 


0.31 


0.50 


yes/yes 


0.19 


0.55 


0.73 


0.25 


0.98 


0.34 


0.26 



KEK 0.23 ± 0.08 0.45 ± 0.10 0.68 ± 0.13 0.27 ± 0.13 0.95 ± 0.04 0.51 ± 0.13 ± 0.05 0.29 ± 0.13 
FINUDA 0.24 ±0.10 



Despite this agreement, we have to admit that more refined and systematic theoretical 
studies should be performed before one can reach definite conclusions from the comparison 
between theory and experiment. For instance, the result obtained for Tnm requires a com- 
ment on the eventual inclusion of the full set of diagrams stemming from the amplitudes in 
Figs. 3 and 4 and eventually from other amplitudes. At first glance, one may think that the 
final outcome from all these diagrams would be a bigger value for r NM , thus spoiling the 
good agreement with data of the present result. This is not necessarily the case, for two rea- 
sons. First, the amplitudes (dl) and (d2) in Fig. 3 and the amplitude (c) in Fig. 4 originate 
from lAlp2h GSC. The inclusion of these correlation amplitudes requires the introduction 
of new terms in the ground state normalization function (5); this leads to a reduction of the 
individual values for each decay width, including the ones we have obtained above. From the 
previous studies in Refs. [18, 19] one observes the following property, introduced by ground 
state normalization: a certain redistribution of the total non-mesonic decay strength among 

24 



the partial contributions occurs when new self-energy terms are included. Secondly, the 
presence of several additional self-energy diagrams which are interference terms between 
amplitudes could also bring to a reduction of the decay rates Tx and T 2 . 

B. Medium and heavy hypernuclei 

In order to have a further indication of the reliability of our framework, which adopts 
the local density approximation to obtain results for finite hypernuclei, we have extended 
the calculation to medium and heavy A hypernuclei. All the GSC contributions and the 
antisymmetrization terms discussed in detail for \ 2 C have been taken into account. The 
results we have obtained are given in Table V and are compared with recent data in Figure 9. 

The GSC-free rate r? represents 86% of the rate r x = T? + r°~ GSC + if sc for \ 2 C. For 
increasing hypernuclear mass number A, this contribution decreases and reaches 81% for 
^ 08 Pb. As expected, GSC contributions are thus more important for heavy hypernuclei. 

The one- and two-nucleon induced rates increase with A and rapidly saturate. Saturation 
is expected to begin for those hypernuclei whose radius becomes sensitively larger than the 
range of the non-mesonic processes. The fact that for fCa and f 8 Ph we obtain very similar 
predictions informs us that in 2 v 08 Pb the non-mesonic decay (both one- and two-nucleon 
stimulated) involve the same nucleon shells which participate in the decay of ^°Ca. Indeed, 
the A wave function (s level of the A-nucleus mean potential) is well overlapped to the 
hypernuclear core already in ^°Ca. 

It should be noted that the slight decrease of the non-mesonic rate r NM going from 
89 Y to \ 39 La is due to the special value of the oscillator parameter fvjj adopted for this 
hypernucleus. Such a parameter, which is obtained as the difference between the measured 
s and p A energy levels in A 39 La, is indeed smaller than the values measured for the two 
neighboring hypernuclei of our calculation, 89 Y and 2 Y 08 Pb. 

The contribution of the two-nucleon induced width is almost independent of the hyper- 
nuclear mass number and oscillates between 22 and 26% of r NM . We note from Figure 9 that 
the datum recently determined at KEK, T 2 = 0.27 ±0.13 [26], is well reproduced by our cal- 
culation. Also the recent determination obtained by FINUDA [27] of F2/Fnm = 0.22 ± 0.08 
for hypernuclei from ^He to \ 6 is in agreement with our predictions. 

Concerning Tnm, the agreement of our predictions with data is also rather good. The 
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only exception is the large underestimation of the datum for the A ~ 200 region, which 
however is also difficult to reconcile with the decay rate measured at KEK for A 6 Fe. No 
known mechanism can be responsible for a large increase in the non-mesonic decay rate 
when going from A 6 Fe to the A ~ 200 region. Concerning the datum for A ~ 200, we have 
to note that, given the difficulty in employing direct timing methods for heavy hypernuclei, 
it has been obtained in experiments (performed at COSY, Juelich [28]) which measured the 
fission fragments (which are supposed to be generated by the non-mesonic decay) emitted by 
hypernuclei produced in proton-nucleus reactions. Large uncertainties affect such delayed 
fission experiments, because of the limited precision of the employed recoil shadow method. 
The produced hypernuclei cannot be unambiguously identified with this method. It is also 
possible that mechanisms other than the non-mesonic decay contributed to hypernuclear 
fission in these experiments. The datum reported in Figure 9 has been obtained as an 
average from measurements for hypernuclei produced in proton-Au, proton-Bi and proton- 
U reactions. 

TABLE V: Decay rates predicted for medium to heavy hypernuclei. 



Hyper nucleus 


L? 


Li 


L 2 


Lnm 


A B 


0.56 


0.64 


0.18 


0.82 


A u 


0.62 


0.73 


0.25 


0.98 


fAl 


0.80 


0.94 


0.28 


1.22 


fSi 


0.81 


0.96 


0.29 


1.25 


A °Ca 


0.87 


1.03 


0.29 


1.33 


| 6 Fe 


0.88 


1.06 


0.33 


1.39 


89 V 
A 1 


0.87 


1.06 


0.33 


1.39 


A 39 La 


0.86 


1.04 


0.32 


1.36 


208p b 


0.86 


1.06 


0.34 


1.40 



We think that the results of the evaluation for medium and heavy hypernuclei are encour- 
aging: they give us some confidence in using the local density approximation for obtaining 
results in finite hypernuclei, even in light systems such as A 2 C. 
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FIG. 9: The predictions for the decay rates T±, T2 and Tnm = ^1 + ^2 are given as a function of 
the hypernuclear mass number A. The results for Tnm are compared with experimental data for 
[29], fC [30], f Al [29], ^ 8 Si [29], fFe [29] and for the region of A between 180 and 220 [28]. 
The datum for T2 is from Ref. [26]. 



C. Closing remarks 

Before concluding, we make here some further comments on our calculation. Through 
our work we wish to emphasize the importance of a detailed many-body treatment of non- 
mesonic decay. This requires the identification and evaluations of a large number of dia- 
grams, working on a step-by-step basis with the perspective of reaching the condition in 
which the terms that are not taken into account can be safely neglected. Considering the 
evolution in the predictions obtained in recent works (see especially Refs. [18, 19]) and here, 
this stability of results has not been achieved yet, and new many-body terms must be con- 
sidered. In our opinion, one should explore the dependencies of predictions on the weak 
transition potential model only after these complicated many-body aspects are properly 
understood. Finally, one should attempt to reach a detailed agreement with experiment 
for Tnm, r„/r p and F 2 /Fnm and thus extract sensible information on strangeness-changing 
baryon interactions. From the experimental side, new and improved data are expected from 
FINUDA@Daphne [31], JPARC [32, 33] and GSI [34]. A direct experimental identification 
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of the two-nucleon induced channels together with the measurement of T 2 is a question of 
particular importance. 

We end this Section with a comment to emphasize the importance of evaluating exchange 
terms. In our many-body inspired calculation, such terms are considered together with GSC 
contributions, which are included on the same ground for one- and two-nucleon induced 
decays. GSC and exchange terms improve by 10% the value of T n /T p . Once GSC are 
included, antisymmetrization turns out to be particularly important for both the one- and 
the two-nucleon induced channels, reducing T 2 by 18% and increasing Ti by 20%. It would 
thus be pointless to neglect exchange terms and evaluate only direct ones. Although the 
introduction of antisymmetry is a difficult task in a many-body framework, one should 
evaluate all those exchange diagrams which are companions of a direct diagram which one 
knows to be relevant. 

VII. CONCLUSIONS 

In this contribution we have studied the effects of GSC in the non-mesonic weak decay of 
A hypernuclei. A non-relativistic nuclear matter scheme has been adopted together with the 
local density approximation, for calculations in hypernuclei ranging from j^B to ^ 08 Pb. All 
isospin channels contributing to one- and two-nucleon induced decays have been considered. 
The employed weak transition potential contains the exchange of mesons of the pseudoscalar 
and vector octets, n, rj, K, p, oj and K* . The residual strong interaction, responsible for 
GSC, has been modeled on a Bonn potential based on ir-, p-, a- and w-exchange. 

By using the Goldstone diagrams technique, GSC have been introduced on the same 
footing for one- and two-nucleon stimulated decays. The normalization of the hypernuclear 
ground state introduced by GSC has been taken into account. We have devoted particular 
attention to those GSC affecting the decay widths T n and T p . The many-body A self-energy 
terms we have considered are originated by the transition amplitudes (a) and (61) of Fig. 3 
(for one-nucleon induced decays) and by the amplitude (a) of Fig. 4 (for two-nucleon induced 
decays). Our approach embodies fermion antisymmetry, i.e., both direct and exchange 
interactions are considered in the various diagrams. Concerning one-nucleon induced decays, 
we have evaluated GSC-free rates generated by amplitude (a), purely GSC terms 

produced by amplitude (61), and interference terms r°"7^ sc between uncorrelated and 
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correlated hypernuclear ground states, i.e., between amplitudes (a) and (61). 

The dominant contribution to I\ = r? + r?" 030 + Tf sc turned out to be r? = T° a + T°. 
For i 2 C, r°- GSC = T n - GSC + r° p - GSC and Tf sc = rf c + rf c represented 13% and 1% of the 
rate r\, respectively; GSC are thus responsible for 14% of the one-nucleon induced width 
(such contribution increases up to 19% for ^ 08 Pb). The above results justify the fact that 
we have neglected the exchange terms in r ^. Exchange contributions are rather relevant 
in the calculation of r°~^j sc (for \ 2 C, they increase r°" GSC by 34% and decreases r°~ GSC by 
8%), while only scarcely contribute to r° ^. GSC and exchange terms together increase the 
value of T n /T p for ^ 2 C by 10%. Thanks to the opening of the two-nucleon induced channel, 
GSC produces a sizable increase (of 32% for \ 2 C when exchange terms are included) in the 
value of T NM = Ti + T 2 . 

The agreement among our final results and recent data is quite good and clearly demon- 
strates the necessity of including GSC and antisymmetrization effects. Nevertheless, we 
believe that a refinement of the present scheme must be pursued. Additional many-body 
terms should be considered, involving for instance the A(1232) resonance. Only after a 
certain stability of predictions is reached within such a microscopic approach one should 
explore the dependencies on the weak transition potential model and determine, through 
detailed comparison with experiment, sensible information on strangeness-changing baryon 
interactions. 
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Appendix A 

In this Appendix we present explicit expressions for the decay rates ^ with PQQ' ^ 
DDD associated to the Goldstone diagrams of Fig. 8 and contributing to Eq. (14). In the 
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main text, these widths have been written as: 

T PQQ' _ \^ q-PQQ' V PQQ' (u u\ (nq) 

L n(p) ~ 'r A ,T A r N ,n(p) L t a ,t a t n K K , K F ) » l Zy J 

T A l,T A ,T N =0,l 

where 

T^ TN (k,k F ) = *\k F )\^(G F rniy^ (30) 

X l l l d V dhdh ' W r A Q t N M^) 

xQ(k,q,q',t,h,h',k F ) - PQQI 5(q - ( £ h'+q~ £ h'^ • 

~ £ 2p2h 

The isospin index ta (t A /) of the weak transition potential is associated to an energy- 
momentum q (q f ), while the nuclear strong interaction isospin index is t n and the correspond- 
ing energy-momentum t. In the following subsections we give the functions W^^ TN (q, q', t) 
and Q(k,q,q' ,t,h,h' ,k F ), the energy denominator an d n (the number of crossing 

between fermionic lines) for the various cases. Finally, we show the isospin sums of Eq. (29). 



i) r D Pf 

' n(p) 

The Wr A ? TATN (q, q', t) function, where q 1 = q and t = h! — h + q, is identical to the 
S^ NT (q,q',t) function in Eq. (A.l) of Ref. [21]. Moreover: 

G(k,q,q',t,h,h',k F ) = 9(q )9(\k - q\ - k F )9(\q - h\ - k F \) (31) 

x6(k F - \h\)6(\q + h'\ - k F \)6{k F - \h'\) , 

£ 2 P 2h = £ 2 P 2h = k o ~ ?k q + £ h q ~~ £ h ~ Vn > ( 32 ) 
and n — 0. The isospin sums are given by: 

r DDE _ -pDDE _| -pDDE , <ypDDE , r DDE _ r DDE , r DDE 

1 n ~~ 1 111 ~r 1 000 T" Ol 101 T 1 110 1 011 t 1 100 

, n-pDDE , -pDDE 
"r oi 001 ~r 1 010 



-pDDE _ kt dde _1_ r DDE _ Q~pDDE , p-pDDS , r DDE _ r DDE 
L p — 01 111 ~r 1 000 01 101 "T 01 110 t 1 011 1 100 

■^DDE _ -pDDj 
001 1 010 



iqrOBE _ -pDDE 
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u ) i „(p) 

The W® A f T ® TN (q, q', t) function, where q' = k — h and t — q, is identical to the S^ NT (q, q', t) 
function in Eq. (A.3) of Ref. [21]. Moreover: 

e(Jfe, q, q', t, h, h' , hp) = 0(q )9(\k - q\ - k F )6(\q - h\ - k F \) (33) 

x0(k F - \h\)6(\q + h'\ - k F \)6(k F - \h'\) , 

F DED _ DDD ( oa\ 
E 2p2h — £ 2p2h ? K^J 

and n — 0. The isospin sums are given by: 

V DED _ _ r DED i V DED , -pDED , n-pDED _ -pDED , -pDED 
L n ~~ 1 111 000 101 oi 110 1 011 "r 1 100 

i opDED , n-pDED 
" t "' 3i 001 "r 01 010 5 

■pDED _ e-pDED , -pDED , rpDED _ opDED , -pDED _ -pDED 
L p — 01 111 "T 1 000 oi 101 oi 110 "T 1 011 1 100 

-pDED _| n-pDED 

1 001 t 01 010 



The VV^^ rjv (?, i) function, where g' = k — q — h! and £ = g, is identical to the 
S^ NT (q,q',t) function in Eq. (A.3) of Ref. [21]. Moreover: 

G(k,q,q',t,h,h',k F ) = 6(q )6(\k - q\ - k F )6(\q - h\ - k F \) (35) 

x6(k F - \h\)0(\q + h'\ - k F \)6(k F - \h'\) , 

EDD _ DDD /o fi \ 

E 2p2h — E 2p2h ; I 00 J 



and n = 1. The isospin sums are given by: 



-pEDD _ _-pEDD i -pEDD _ -pEDD , typEDD , -pEDD , <rpEDD 
1 n — 1 111 ~r 1 000 1 101 61 110 "r 1 011 61 100 

, pECB , -pEDD 
' 001 010 ' 

-pEDD _ _r^-pEDD , -pEDD , -pEDD _ typEDD , rpEDE , opEDi 
1 p — 01 111 1 ooo ~r 1 101 01 110 ~r 01 011 01 100 

-pEDD _ -pEDD 

1 001 1 010 
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iv) r™ 

The W^f T ^ TN (q,q',t) function, where q' = k — h and t = h — hi — q, is identical to the 
S^ NT (q,q',t) function in Eq. (A.7) of Ref. [21]. Moreover: 

e(Jfc, q, q', t, h, ti, hp) = 0(q )6(\k - q\ - k F )9(\q - h\ - k F \) (37) 

x6{k F - \h\)6(\q + h'\ - k F \)6(k F - \h'\) , 

DEE _ DDD /oo\ 
E 2p2h — e 2p2h J \ M ) 

and n — 1. The isospin sums are given by: 

■pDEE _ rpDEB _| pDEE , pDEE , pDEE , rpDEB , pDEE 
1 n OX 111 ' 000 101 "r 1 110 di 011 1 100 

,-pDEE | -pDEE 
~T~ L 001 "T 1 010 ) 

P_D_B_B _ rf-pDEE a pDEE , < r DEE , r,-pDEE , pDEE 
L p — ZL 111 ^ 101 ^ L 110 "r Zi 011 Zi 100 

"r zi 001 "r Zi 010 



\ r EDE 
V J 1 n(p) 

The W^f^ Tjv (g, g', t) function, where q' = k — q — hi and t = hi — h + q, is identical to 
the Sp? NT (q,q',t) function in Eq. (A.7) of Ref. [21]. Moreover: 

G(k,q,q',t,h,h',k F ) = 6(q )6(\k - q\ - k F )6(\q - h\ - k F \) (39) 

x9(k F - \h\)0(\q + h'\ - k F \)6(k F - \h'\) , 

EDE _ DDD 

E 2p2h — E 2p2h ; l^ U J 

and n — 1. The isospin sums are given by: 

pEDE _ _-pEDE i -pEDE , opEDE , pEDE _ pEDE , pEDE 
1 n — 1 111 "T 1 000 "r 01 101 ~r 1 110 1 011 1 100 

, npEDE , pEDE 
~l~ oi 001 1 010 

pEDE _ ApEDE I apEDE A pEDE cypEDE , rypEDE , rypEDE 
1 p _ '* i 111 "T 01 101 ^ 110 Zi 011 "T Zi 100 "T Zi 010 
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vi) T™» 

The W^f r ^ Tjv (q, q', t) function, where q' = k — h and t — k — q — h! , is identical to the 
S^ NT (q,q',t) function in Eq. (A.5) of Ref. [21]. Moreover: 

Q(k, q, q', t, h, ti, k F ) = e(q )9(\k - q\ - k F )9(\q + h + ti - k\ - k F \) (41) 

x0(k F - \h\)0(\q + h'\ - k F \)6(k F - \h'\) , 

4v2h =k -e h + e q+h+h > k - e q+h> - V N , (42) 
and n — 1. The isospin sum are given by: 

-pEED _ _-pEED i -pEED , -pEED , orEBD _ -pEED , -pEED 
1 n ~ 1 111 "r 1 000 101 oi 110 1 011 ~r 1 100 ' 

i pEED | o-pEED 
' 001 "T 01 010 

-pEED _ A -pEED A-pEED < (ypEED n-pEED , n-pEED , n-pEED 

1 p — in ioi -r oi no zi on -I- zi 100 -|- zi 001 . 



vii) r EEE 

7 n(p 



The W^f T ^ riv (g, g', t) function, where q' = k — h and t = h + q — k, is identical to the 
S*?* (q,q',t) function in Eq. (A. 9) of Ref. [21]. Moreover: 

Q(k,q,q',t,h,h',k F ) = e(q )e(\k-q\-k F )e(\h' + q\-k F \) (43) 

x6{k F - \h'\)9{k F -\k-h + h'\)6{k F - \h\) , 

=k -e h + e h - e k h+K - V N , (44) 
and n — 0. The isospin sums are given by: 

-pEEE _ trpEEE , -pEEE , -pEEE , -pEEE , tr-pEEE , -pEEE 
1 n — 01 111 ' 000 101 1 110 ~r di 011 ~r 1 100 

i^t^E E E t^E E E 
~* 001 "T 1 010 ) 

EEE wEEE , tEEE , rr^EEE , r F EEE , -rEEE tEEE 



i p — a 111 + i 000 + 01 101 ~r 01 110 + i 011 ~~ 1 



100 



t^E E E t^E E E 
~ L 001 — 1 010 
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Appendix B 



The explicit expressions of the function TZ(q , q) of Eq. (22) reads: 



K(q ,q) = 



7T 



m 



(2tt) 3 q 



x In 



2mq 



2qkj 



r 



q 



2m j m 



k F + 





2mq — q 2 + 2qk F 



+ d(2k F -q) 



k> 



-f) : 



In 



2mq + q 2 — 2qk F 



2mq — q 2 + 2qk F 



- '\ (2 9o 3. (*, - i) 
q l \ m V 2/ 

-k F ) 



x In 



2mg — q 2 + 2qkp 



where q 



2mq + q 2 — 2qk F 
\q\ and m is the nucleon mass. 



q Q m In 



go«i 2 - q 2 {k F - q/2f 



m 2 q^ 
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